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The partition function for the free neutral scalar field given in the framework of the method of the 
path integral quantization was exactly derived by the new algebraic method. We found that the free 
massive neutral scalar field on the hnite lattice is thermodynamically inconsistent in the whole range 
of the variables of state T and V including the thermodynamic limit because its thermodynamic 
potential is not a homogeneous function of the hrst order with respect to the extensive variable of 
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consistent either at high temperatures or in the thermodynamic limit. At Hnite values of volume 
the free neutral scalar Held in the continuum limit is thermodynamically inconsistent. 
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I. INTRODUCTION 

The process of particle (hadron) production from the 
hot medium created in heavy ion collisions and hadron- 
hadron reactions at the LHC dHi and the RHIC 
calls for the introduction of finite temperature and vol¬ 
ume into studies of quantum field theory (QFT). How¬ 
ever, most theoretical approaches are used to describe 
the thermodynamic properties or/and phase transitions 
for relativistic quantum fields (i.e. bosons and fermions, 
which are generated in these processes) are considered 
in the thermodynamic limit. In these models, the infi¬ 
nite volume of the system is implemented directly 
or indirectly through the application of the continuous 
momentum spectra instead of the discrete eigenvalues of 
the momentum operator of particles [iM3- Thus, the 
finite size effects are important for consideration in the 
statistical theories meant to describe the nucleus-nucleus 
and proton-proton collisions. They may significantly in¬ 
fluence the thermodynamic characteristics of quantum 
fields and, in particular, of the QCD phase diagram as a 
whole [l^ . 

There are various methods of quantization of fields, 
such as the canonical method, the path-integral method, 
etc. (Tgl - f^ . The thermodynamic properties of the QCD 
phase diagram are described by both the canonical quan¬ 
tization formalism (second quantization) [T^ . and the 
path-integral formalism [23|. Numerical computations 
in the framework of the path-integral formalism by the 
Monte Carlo method, for example in the case of lattice 
gauge theory, allow the study of nonperturbative effects, 
for example those manifesting themselves in the QCD 
phase transitions . All of these methods of field quan¬ 
tization are equivalent BUI. 

The formalism of statistical mechanics agrees with the 
requirements of equilibrium thermodynamics, i.e. it is 
thermodynamically self-consistent under the condition 
that the thermodynamic potential of the statistical en¬ 
semble, which contains all information about the physical 
system, is a homogeneous function of the first order with 
respect to the extensive variables of state [2^ - 1^ . For 


most physical systems this condition is fulfilled in the 
thermodynamic limit. The problem of thermodynamic 
consistency in the case of lattice quantum chromody¬ 
namics (LQCD), which is hoped to accurately describe 
the phase transitions for the relativistic quantum fields 
featuring in heavy ion collisions, has not been rigorously 
studied yet. It is considered that the trace anomaly en¬ 
hancement is the interaction measure in both the pure 

S age theory on the lattice [3l| and the full LQCD [3^ 
. In the approximate numerical calculations for the 
lattice QCD as given, for example, in [33l - l35j| it is difficult 
to verify the homogeneity properties of the thermody¬ 
namic potential of the system. Therefore, to investigate 
the severity of this problem, an exactly solvable statis¬ 
tical model for a free relativistic quantum field in the 
framework of path-integral quantization on a finite lat¬ 
tice should be reconsidered and finite volume effects as 
well as the continuum limit should be carefully analyzed. 

The main purpose of this paper is to verify the thermo¬ 
dynamic self-consistency of the free neutral scalar field on 
the lattice and in the continuum limit in a finite volume. 

The structure of the paper is as follows. In Sec¬ 
tion El we briefly describe the path-integral formalism 
for the neutral scalar field in the configuration and mo¬ 
mentum spaces. The continuum limit for the thermody¬ 
namic quantities is given in Section IHII The vacuum and 
physical thermodynamic quantities on the finite lattice 
are defined in Section IlYl The condition of the thermo¬ 
dynamic self-consistency of the statistical model in the 
grand canonical ensemble is formulated in Section|Vl The 
results are discussed in Section IVll The main conclusions 
are summarized in the final section. 


II. PATH INTEGRAL QUANTIZATION FOR 
FREE REAL SCALAR FIELD 

The method of path integral quantization for the real 
scalar field is well known. We propose a new algebraic 
method to solve the one-spatial dimensional partition 
function and recite the most important path integral re- 
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lations to keep the paper self-contained. We consider a 
system consisting only of the neutral scalar field (/>(a;) 
with mass m. The corresponding classical action S is the 
integral over space-time 

5” = y dtL = J d'^xC{<p{x),d^<p(x)), (1) 

with the classical Lagrangian density C given by 

^ ( 2 ) 

The thermodynamics of the real scalar field of volume 
V, in contact with the heat and particle reservoir of tem¬ 
perature T and chemical potential fi, is defined by the 
statistical operator and the partition function, 

g = and Z = Tr (3) 

Zi 

respectively, where /3 = 1/T, H is the Hamiltonian of 
the system and Q is the electric charge operator in the 
canonical operator formalism. For the real scalar field 
the charge operator <5 = 0- 

The partition function in the canonical operator for¬ 
malism may be redefined in the form of the functional 
integral with the classical action. This method gives an¬ 
other consistent description of quantum field theory . 

Thus, the partition function ([3]) in the path integral rep¬ 

resentation can be rewritten as 

Z = (4) 

The field operator (l>{x) in the Schrodinger picture does 
not depend on time t, 4){x)^ and has the eigenstate \4>{x)) 
with the eigenvalue </>(a?): 

= (j){x)\(j){x)). (5) 

Similar relations are satisfied by the canonical conjugate 
field operator tt{x) = dC/d(j){x) = (l){x). The Hamilto¬ 
nian H of the system in the Schrodinger picture is ex¬ 
pressed in terms of the field operators ^(x) and 7r(x) as 


and volume V form a 4-dimensional lattice A, the cells of 
which are given by the vector = (n, np). Then to each 
cell there is attributed the field operator 4>{n^) with 
an eigenstate |(/)(n^)) and eigenvalue (j){n^). We have a 
discrete set of a finite number of operators. For con¬ 
venience, we denote the field operator in another form 
= (j){n,j) {I € A 3 ,i = np), where I is an integer and 
A 3 is the 3-dimensional subspace of the lattice A. 

For simplicity, to start with, let us find the analytical 
expression for the partition function in one spatial di¬ 
mension and generalize it finally to the usual three spa¬ 
tial dimensions. The partition function Q in one spatial 
dimension can be written as 

Z = lim Ziat, ( 8 ) 

N(j—¥(yD,N^—¥cc> 

Zlat= / 

1=1 

^ ^ ( 9 ) 

where lim denotes the continuum limit which means that 
Ncr —?> oOjOct —>■ 0 at Lct = ttcrNcr = coust and Njs —>■ 
oo,a,g—^0at/3 = apNp = const. The exponents in 
Eq. do]) are enumerated from right to left by the index 
i = 1,..., Np. The state vectors \4>i) = |0qi,..., 
and Itt^) = ..., TTN^,i) are orthogonal and complete: 

(jpj), (10) 

1=1 



To solve Eq. ([9]), we insert in the left side of each z-th 
exponent the product of two unit operators for the fields 
\(j)i+i) and \'Ki). Using Eq. (ITUl) . we obtain 


i? = y d^x'H{TT{x),p{x)), 

(6) 

73 = i - 1 - (V(^)^ - 1 - . 

(7) 


The main point of the path integral quantization 
method is that the “time” interval of continually variable 
length /3 is subdivided into Np equal intervals of length 
ap, P = apNp, and the volume V of the system is divided 
into small cells, each with volume so that V = L^., 
where = a^Nu. The small intervals on the /3-axis 
are labeled by the integer np = 1,..., Np and the space 
cells in V are fixed by the integer vector n = (ux, ny^Uz) 
with the coordinates = 1,. ■., ■ The discretized /3 


Np 


z«= /nn 


1^1 2=1 


(XfjdiTI id(l)i i 
27 T 


iVa 


Y\_^4'l,Np + lS{Pl,l — 4'l,Np + l) 


1=1 

N, 


P 

i=l 

The matrix element is 


-a^ap y^XT-Lip-pQip] 

= e ‘=1 (TTiki), (15) 
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where the Hamiltonian function © for the Klein-Gordon 
field on the lattice in one spatial dimension is given by 


+ 2 


2 I 1 <Pl,i \ -L 2/2 


+ 2 ”^ fiv 


(16) 


Note that for the real scalar field the chemical potential 
/r = 0. 

Using the eigenstate Itt^) in the field representation, 
the Hamiltonian function for the Klein-Gordon field on 
the lattice rewritten for one spatial dimension (HU) 
and the (anti)periodic boundary conditions for the field 
in the (single) spatial direction, we can write 


Zlnt — 


ITIT 

1=1 i=l 


27r 


1=1 

Np 

d(j)N, + l,i5{4’l,^ — i<y(l>N„ + l,i) 


exp 


N„ Np 

--EE 

1=1 i=l 


,( 17 ) 


where Co- = 1 and Co = — 1 for the periodic and antiperi- 
odic boundary conditions along the spatial x-axis, respec¬ 
tively. Integrating over the variables 7 r/^i, (jii^Np+i and 
4'N,+i,i under the condition that Re(ao-a/ 3 / 2 ) > 0 and in¬ 
troducing the new variable (j)j{i,i) = (ao/ 2 a/ 3 )^/^^iy with 
J{1, i) = Np{l — 1) -b z, we obtain 


N„ Np 


^lat — 




exp 


—00 


= 1 


-EEEE Aj{lp)J{k,])4>J(Ui)4>J(k,j) 

1—1 i—1 k—1 j—1 


(18) 


and 


Ajj' = Si^kSij 




dkk [(1 “t“ (1 dj^Np') 

— [di,lSj^Np + di^NpSj^l] 

~ di,j [(1 — 5k,1+1 + (1 — 5k,N^) 5l,k+l] 


— 5i,j ^<j{5i,i5k,N^ + 5i,N^5k,i) 


(19) 


where Ajji = Aj(i^i'jjf^k,j)- The matrix H is a symmetric 
square matrix. Thus, integrating (fTSl) with respect to the 
variables 4’J(i,i ), obtain 


^lat — 


■\/det A 


( 20 ) 


The determinant of the matrix A can be calculated by 
the recurrence equations of the form 


NxNp 


detA= pi, 

1^1 

( 21 ) 

Pi = Au-Y.^^ 

( 22 ) 

2—1 

n --9A - V 

q., - ZA,, ’ 

(23) 


where Re(p;) > 0, ^ = 1,..., N^Np, i = 1,..., N^Np — 1 
and j = z -b 1,..., NaNp. The matrix elements Aij of the 
matrix A are given by Eq. (HU). 

Let us rewrite the partition function ([^ in the momen¬ 
tum space on the basis of the Fourier transform defined 
on the lattice. The four-dimensional momentum space 
A, which corresponds to the lattice A, can be defined as 

~ Stt 

A = {P/. = {Px,Py,Pz:Pp)\Pp. = -+ dp), 

kpi<kp<kfj,2,P = x,y,z,l3} (24) 


and 


kpi = - 


1 I Vp^p 


kp2 = 


2 2 ’ 
Np-^ , Vp^p 


(25) 

(26) 


2 2 ’ 

where C/x = 1 (d/i = 0 ) and C^t = “1 {dp = 1 / 2 ) for the 
choice of periodic and antiperiodic boundary conditions, 
respectively, along the /z-direction and 77 ^ = 1 for 
even and rjp = 0 for odd [^ . 

It is not difficult to prove that the vectors and 
satisfy the following relations: 




dk h' — - 

_/y 




(27) 


N,, 


Tin, — ! 


E 


a /1 (28) 


and 


1 

_ = A, 

N,, 




riu^l 


[(1 20^)77^-b 20^(1 77^)]) (29) 

where -{N^ - 1) < k^ - k'^ < - I and -{N^ - 1 ) < 

Ufj—n'^ < fV^ —1. Then the Fourier transform for the field 
4>J = 4'J{na,,np) on the lattice with one spatial dimension 
and its inverse transform can be written as 

kx2 ^02 r 

(k = pi{aaPxnx+a^P0n^) 

kh..kp=kp,V^^p 


N„ Np 


/^=E E 


0J 




^ — i{aapxnx+apppnp) 


(31) 


rix — l 
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where the index J = J{nx,np) and the index / = 
I{kx,kij) = Njjikx — kxi) + kp — kpi + 1. The neutral 
scalar field is real and (jfj = (j)j. Therefore, the complex 
function // = fi^k^ykis) can be represented by its ampli¬ 
tude and phase in the form 

(32) 


Substituting Eqs. (1321) and (l30l) into and using 
Eq. (l29ll . we can write 


7 ^ 

Ziat = '^ " IdetJ^I / dRi 

_00 fc/3=fef)l 


exp 


kx2 ^/32 kx2 

- E E E E Bn,RjRj, 

kx=kxl fe/3=fe/31 k'^=kxl fc^=fc/31 


(33) 


and 


Bii' = Gi [Skp+k'^,0 + Vd^kfi+k'^,Np] 

\^kx-\-k'^-\-‘ 29 x ^0 T C , 7 Vct ] J 

c = (l-20,)ry, +20,(1-r/,), 


Jji = 


1 


d(l>j _ _ 

dRi ~ 


c,i[axPx{nx + l)+ai}pp(np+l)\ 


Gj=aj 


2 . appp \' 

— sm 

0/3 2 / 


2 . (kfjPx \ 

— sm —-— 

Oo- 2 / 


I 2 2 

“h QaTTl . 


(34) 

(35) 

,(36) 


(37) 


The matrix _B is a symmetric square matrix. Thus, inte¬ 
grating (1^ with respect to the variables Rj, we obtain 


_ I detj"! 

Vdet B 


(38) 


The matrix B has Na-Np nonzero elements equal to Gj. 
In any row and any column there is only one nonzero 
element, regardless of the periodicity conditions for the 
cc-axis and parity conditions for and Np. Therefore, 
the determinant of the matrix B can be written as 

kx2 ^92 

detB = x n n Gi, 

kx—kxi kjs—kjsi 

where Ip = {I - pp)3Np + pp{Np + 2) and /, = (!- 
77 ,) 20 ,-|- 77,(1 —20,). Substituting Eq. (1^ into Eq. ([38|) . 
we obtain 


(39) 

(40) 


Ziat = . E ^=, (41) 

Y\kx2 TT«/32 ^ 

where \ dei J\x~^^'^ = 1 and Gj is the function calcu¬ 
lated by formula (I37p . Note that Eq. (HD) is equivalent 
to Eq. (1^ . 


Now the partition function ()41l) for one spatial dimen¬ 
sion can be generalized to three spatial dimensions of the 
momentum space. For this reason, we can rewrite the 
partition function m and the function (1371) in the form 



where 9p = 0, p /3 = 2Trkp/{apNp) and a = 1,2,3. Let 
us remark that the lattice partition function for the 
particular case of the periodic spatial boundary condi- 
tions, 0 O- = 0 , was firstly obtained by another method 
in [ 11 ]. 

The thermodynamic quantities on the lattice are de¬ 
rived from the partition function (1421) . The density of 
the thermodynamic potential wg = —{f5V)~^ hiZiat-, the 
energy density Ee = —V~^d\nZiat/dl3, the pressure 
Pe = Ziat/dV, and the entropy density se can 

be written as 


Ue 


£e 


Pe 


fe/32 


2PV 




k=k^i kp=kpi 

kx2 ki32 2 

OJ 


1 




k—kfji ^/3i 

. k^2 ^02 2 2 

1 CJ — m 


E E 


3/3V ^ ujj+uj'^ 


(45) 

(46) 

(47) 


and 


Se = fi{—i^E + £e), (48) 


where dNp/d(3 = 0 and dN^/dV = 0. Note that the en¬ 
ergy density (l46l) and pressure (iTTl) with the periodic spa- 
tial boundary conditions, Oa = 0, were obtained in |36l |. 


III. CONTINUUM LIMIT 

Let us find the density of the thermodynamic poten¬ 
tial (l45l) . the energy density (l46l) and the pressure (iTTl) 
in the continuum limit. Here we consider only the pe¬ 
riodic boundary conditions, 0cr = 0 and 0p = 0, as in 
Ref. [ 3 ^. In this case, we have Pa = 2'nka/{aaNa) and 
the functions (l43l) and (l44l) can be written as 


2Np 

. Ekp 

Np ’ 

(49) 


'2Na- . TTkaY , 2 

—— sm -1- m , 

^ La Na J 

(50) 
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where = —(7V^/2) + 1, k ^2 = N^/2 for even and 
k^i = -{N^ - l)/2, k ^2 = {Nfi - l)/2 for odd, and 
a = 1,2,3. 

First, let us find the thermodynamic quantities (HSl) - 
m in the limit N/s —>■ oo, O/j —>■ 0 at /3 = apNp = const 
and fixed value of N^. Expanding the function under the 
sums in Eq. (1461) into the series on around the point 
Np = 00 and taking the limit Np —> oo, we obtain 


-I k^2 oo r o \2 

J_ V V (/^^) 


(51) 


Erom now on the asterisk symbol will denote the thermo¬ 
dynamic quantities given in the limit Np —>■ oo, ap —)> 0 
at /3 = const and fixed value of N^. Using the generating 

function (see, for example, [13) 


oo 9 

E a^ a a 

{2'KkY + a? ^ 2^°^^r 

k— — oo 


(52) 


one obtains 


, 1 Puj Buj 

= ^ ^COth— =£„ + i 

k=k„i 




k—k„ 

ktT2 




V - 1 ’ 

k-kn-l 


(53) 

(54) 

(55) 


where e* and e* are the vacuum and physical terms, re¬ 
spectively, of the energy density e*^. For m = 0, the zero¬ 
mode term k = 0,kp = 0 under the sums in Eqs. (H51) . 
dSU and the zero-mode term k = 0 under the sums in 
Eqs. (15^ . (1551) are suppressed 
The pressure (H71) in the limit Np —>• oo and constant 
value of AU can be written as 



where p* and p* are the vacuum and physical terms, re¬ 
spectively, of the pressure p’^. Note that for m = 0 the 
zero-mode term k = 0,kp = 0 under the sum in Eq. (iTfl) 
and the zero-mode term fc = 0 in Eqs. dSTl) and (1551) are 
all equal to zero [s^. Thus, for the free massless real 
scalar field we have the relations pe = £e/^, Pe = 
pI = e*/3 and p* = e*/3. 

Considering Eqs. (1551) . (|55)) and the definitions of the 
quantities (|55 1) -(|T7 1) . the density of thermodynamic po¬ 
tential (|551) in the limit Np —> oo and fixed value of AU 


can be written as 


a;|=w:+w;^. 


= 


, k„2 

V E 


2U 

1 


w, 


k — fcjy 1 
k(j2 


= ^ Y. In (1 - e , 


(59) 

(60) 

(61) 


k—kfj 


where w* and are the vacuum and physical terms, 
respectively, of the density of the thermodynamic po¬ 
tential Note that the functions (051), (E® and dSD 
in the case of the massive neutral scalar field converge 
at Np —>■ oo and constant N„. The numerical proof of 
Eq. (1551) will be given in Section |Vll For m = 0, the 
zero-mode term fc = 0, = 0 in Eq. (USD and the zero¬ 

mode term fc = 0 in Eq. m are infinite. These terms 
describe a condensate of massless bosons which should 
be excluded from the thermodynamic potentials because 
this condensate is an unobservable. 

Now let us find the thermodynamic quantities (1451) - 
(|i71) in the continuum limit Np —>■ oo, A^o- —t oo, ap —>■ 0, 
tto- —!> 0 at /3 = apNp = const and = OaN^ = const. 
Then, the energy density in the continuum limit can be 
rewritten as l20l l38ll 


£e — ^ I 


= ^E“. 


c 1 1 

fc 


(62) 

(63) 

(64) 


where oj is the one-particle energy (1501) . It takes the form 


27r 


= -I- m2, p^ = — ka, 

J^rr 


(65) 


where ka = 0, ±1, ±2,.. .±oo and a = 1, 2, 3. For m = 0, 
the dispersion relation (j65|) is a; = |^ and the zero-mode 
term k = 0 under the sum in Eq. dSl is suppressed as 
in [ 3 ^ . From now on the index ’c’ will denote the physical 
quantities given in the continuum limit. The pre ssure 
dnoi in the continuum limit can be rewritten as |20l . 1^ 



where oj is given in Eq. (1051) . Note that for m = 0 the 
zero-mode term fc = 0 in Eqs. dlTl), dSll) is equal to 
zero [ 33 . Thus, for the free massless real scalar field 
in the continuum limit we have the relations = £s/3, 
pI = e%/Z and p" = e73. 
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The density of the thermodynamic potential (|5^ -" (l6T1) 
in the continuum limit is given by 


- ^v+ 

(69) 

k 

(70) 


(71) 


k 


where uj is defined in Eq. (I65|) . For m = 0, the zero¬ 
mode term fc = 0 in Eq. (EU is infinite and should be 
suppressed because it is nonphysical. Thus, in the con¬ 
tinuum limit the energy density (I5^ - (|M)) . the pressure 
(l66l) - ((68l) . and the density of the thermodynamic poten¬ 
tial (IS^ - dTTl) exactly coincide with their corresponding 
quantities obtained by the method of second quantiza¬ 
tion [ 2 O, Hence, it was proved that for the free neu¬ 
tral scalar field the method of path integral quantization 
and the method of canonical quantization are equivalent 
at any values of T and V. 

Let us investigate the trace anomaly and the potential 
inhomogeneity for the free neutral scalar field in the con¬ 
tinuum limit. Using Eqs. (I64[) and (1681) we can write the 
trace anomaly in the continuum limit as 

= (72) 

k 

The trace anomaly A 2 in the continuum limit defines 
the deviation of the equation of state from the Stefan- 
Boltzmann equation. The trace anomaly dZlD for the 
massive neutral scalar field in a finite volume and at low 
temperatures is nonvanishing; however, it decreases with 
T. For TO = 0, the trace anomaly dZl is equal to zero, 
A 2 = 0, because the zero-mode term A: = 0 in Eq. dZH) 
was suppressed. The potential inhomogeneity in the con¬ 
tinuum limit can be defined as 

A^ ^ ^ ^ In (1 - e-^^) 

k 



k 


This quantity is a measure of violation of the homo¬ 
geneous properties of the thermodynamic potential of 
the grand canonical ensemble. The system in the grand 
canonical ensemble is thermodynamically inconsistent if 
the quantity (1751) is not equal to zero. For the proof of 
this statement see Section Ivl 

In the thermodynamic limit as U —> 00 , the ther¬ 
modynamic quantities of the massless neutral scalar 
field recover their Stefan-Boltzmann limit, i.e., ess = 
7r2/(30/?^), psB = £sb/S and lusb = -PSB- However, 
at finite values of volume V and temperature T the shift 
Asb = £sb — ioT the massless neutral scalar field is 
not equal to zero and the energy density deviates from 


its Stefan-Boltzmann limit. Note that the trace anomaly 
dZH) for the massless neutral scalar field in the Stefan- 
Boltzmann limit is zero. 


IV. PHYSICAL THERMODYNAMIC 
QUANTITIES ON A EINITE LATTICE 

We define the vacuum and physical terms of thermo¬ 
dynamic quantities on the finite lattice {Np,N^) as in 
Ref. [ 3 ^. In this case, the vacuum terms of the ther¬ 
modynamic quantities on the finite lattice are derived 
from Eqs. dm-dizi) by changing the summation over kp 
to integration, which corresponds to T = 0. Thus, by 
definition, we have the vacuum terms as [ 3 ^ 



where uj is given in Eq. (pl). Then, the physical terms 
of the thermodynamic quantities on the finite lattice can 
be written as [3g 



These quantities are the same as those from Ref. [ 3 ^ . 
Using Eqs. dzsl) and dZH), we obtain the trace anomaly 
for the free real scalar field on the finite lattice as 

k,k(3 ^ 
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It defines the deviation of the equation of state from the 
Stefan- Boltzmann equation, psB = ess/S. The trace 
anomaly for the free massive real scalar field on the finite 
lattice is not equal to zero at the finite temperature and 
volume. However, for m = 0 the trace anomaly (j80|) is 
equal to zero, A2 = 0. The potential inhomogeneity for 
the free real scalar field on the finite lattice can be defined 
as 


Ai = + P) = ^ ^ In [ajiujj + w^)] 

k,k^ 



It defines the violation of the homogeneous properties 
of the thermodynamic potential for the finite values of 
{Np,Ncr). See Section M 

It should be stressed that the vacuum and physical 
quantities ujph, Sy, e and Pv, p on the finite lat¬ 
tice given in Eqs. (TTHl -dT ^ resemble the thermodynamic 
quantities ([70ll . (TtTII . (l63l) . (|64ll and (|67|l . (IMI) . respec¬ 
tively, in the continuum limit. Thus, in the continuum 
limit the vacuum and physical terms of the thermody¬ 
namic quantities given in Ref. exactly coincide with 
the same quantities obtained in the canonical quantiza¬ 
tion method. 


V. THERMODYNAMIC POTENTIAL AND 
ZEROTH LAW OF THERMODYNAMICS 

In the equilibrium statistical mechanics all thermo¬ 
dynamic quantities belong to the class of homogeneous 
functions of the zero and first order. This property of the 
thermodynamic quantities provides the fulfilment of the 
requirements of the equilibrium thermodynamics. For ex¬ 
ample, in the grand canonical ensemble the zeroth law of 
thermodynamics is satisfied if the temperature is inten¬ 
sive (zero order) and the grand thermodynamic potential 
is extensive, i.e. it is a homogeneous function of the first 
order with respect to the extensive variable of state V. 

Let us consider the case when the thermodynamic po¬ 
tential of the grand canonical ensemble is indeed a ho¬ 
mogeneous function of the first order with respect to the 
extensive variable of state V. We have (^.[sOjl 

n{T,V,p) = Vuj(T,fi). (82) 

Then the pressure p and the potential inhomogeneity 


Ai = 3(w +p)/T^ can be written as 

dVL fl 
^ ~V 

and 


Ai 


3 


ji4 



= 0 . 


(83) 


(84) 


For the free neutral scalar field p, = 0. Note that only in 
this section uj denotes the density of the thermodynamic 
potential of the grand canonical ensemble. 

To prove the zeroth law of thermodynamics, let us di¬ 
vide the system into two subsystems (1 and 2). Then the 
extensive variables of state of the grand canonical ensem¬ 
ble should be additive and the intensive variables of state 
should be the same 

V = V^ + V^, T = T^ = T^, p = p} = p'^. (85) 

Since the function uj in Fq. (15^ depends only on the 
intensive variables of state T and p then we have 

u:{T,p)=uj\T\p^)=uj\T\p'^). ( 86 ) 

Using Fqs. (15^ . (1551) and (1551) . we obtain 

n^{T\v\p^)+n^{T‘^,V‘^,p^) = ^{T,V,p). (87) 


Thus, we have obtained that if the thermodynamic po¬ 
tential is a homogeneous function of the first order and 
the temperature is intensive, then the thermodynamic 
potential is an additive function and the potential inho¬ 
mogeneity is zero. This proves the zeroth law of thermo¬ 
dynamics for the grand canonical ensemble [2^ . 

Let us consider a more general case when the thermo¬ 
dynamic potential of the grand canonical ensemble is an 
inhomogeneous function. For example, we can write 


n{T,V,p) = V^Cj{T,p), (88) 


where a is a real number. For a = 1 we have Fq. (1821) . 
Then the pressure p and the potential inhomogeneity Ai 
can be written as 


p = - 


dil 

W 



(89) 


and 


Ai 


3 


y4 



2.4(1 a )^. 


(90) 


For the inhomogeneous thermodynamic potential (1881) at 
a ^ 1 the potential inhomogeneity Ai 7^ 0 and the pres¬ 
sure p ^ —U/U. 

To verify the zeroth law of thermodynamics for the 
inhomogeneous thermodynamic potential (I88L let us di¬ 
vide the system into two subsystems (1 and 2) and require 
Fq. (155)1 . Then we have 

Cj{T,p)=Cj\T\p^)=Cj\T\p'^). 


(91) 













Using Eqs. (ESI), EH) and EB), we obtain 


(]/!)« + (y2)a 


n{T,v,fi). 


(92) 


For a ^ 1 we have U. Thus, we have ob¬ 

tained that if the thermodynamic potential is an inho¬ 
mogeneous function and the temperature is intensive, 
then the thermodynamic potential is nonadditive and 
the potential inhomogeneity is not equal to zero. This 
means that the inhomogeneous thermodynamic poten¬ 
tial of the grand canonical ensemble and corresponding 
relation p ^ —U/U (Ai ^ 0) violate the zeroth law of 
thermodynamics and thus the formalism of this statisti¬ 
cal mechanics is thermodynamically inconsistent. There¬ 
fore, the grand canonical ensemble is thermodynamically 
inconsistent if the density of the thermodynamic poten¬ 
tial is not equal to minus pressure and the potential in¬ 
homogeneity is not equal to zero. 


VI. ANALYSIS AND RESULTS 



10 ' 10 ' 10 ' 10 " 10 ' 


Let us study the thermodynamic properties of the real 
scalar field on the lattice and in the continuum limit at 
a finite temperature T in some finite volume Y which 
are characteristic of the ultrarelativistic heavy ion and 
hadron-hadron collisions 

First of all, let us numerically verify Fq. (EB saying 
that the density of the thermodynamic potential (I45p in 
the limit iVg —>■ oo and fixed value of is divided into 
the sum of two terms: the vacuum term and the physical 
term. Figure [1] represents the density of the thermody¬ 
namic potential wg (liSll . its limit (1551) . and the shift 
= oj'^ — loe functions of A(a for the free neutral 
scalar field on the finite lattice with the mass of 7r° pion 
at the temperature T in some volume V for different val¬ 
ues of N^. For the given T, V and N^, the density of 
the thermodynamic potential lue (HSI) increases with A(a 
and attains the constant uj'^ (1551) at large values of Np. 
The shift A,^ decreases with increasing V/3 and tends to 
zero, A,^ = 0, as > oo, which numerically proves 

that Fq. (1591) contains only two terms: the vacuum term 
w* and the physical term The quantities uje and 

increase with No-. Thus, in the continuum limit, the 
density of the thermodynamic potential oje and tend 
to infinity, but the shift A^ is equal to zero. In the con¬ 
tinuum limit Fq. (159|) recovers Eq. (1551) . Thus, contrary, 
for example, to Ref. po| . the density of thermodynamic 
potential (155)1 for the free neutral scalar field in the con¬ 
tinuum limit does not contain any infinite constant. The 
thermodynamic potential (1551) is composed of the vacuum 
and physical terms only. 

Let us investigate the thermodynamic self-consistency 
of the free massive neutral scalar field on the finite lat¬ 
tice. The formalism of statistical mechanics in the grand 
canonical ensemble agrees with the requirements of equi¬ 
librium thermodynamics if the thermodynamic potential 


FIG. 1. (Color online) The density of the thermodynamic 
potential ljJe (|45]), its limit oj% (|59l). and the shift A„ as func¬ 
tions of Njs for the free neutral scalar field on the finite lattice 
with the mass of 7r° pion at the temperature T = 100 MeV in 
the volume U = 9® fm® for different values of Na. The lines 
1,2 and 3 correspond to the values of No- = 10,20 and 30, 
respectively. 


of the grand canonical ensemble is a homogeneous func¬ 
tion of the first order with respect to the extensive vari¬ 
able of state, the volume V. This means that the density 
of the physical thermodynamic potential uiph (l77l) and 
the physical pressure p (1751) should fulfill the equation, 
ujph = —p, for any values of the variables of state T and 
V, i.e., the potential inhomogeneity Ai (I5T1) should be 
equal to zero [2Bj^ . See also Section [V) 

Figure [2] represents the behavior of the physical en¬ 
ergy density dZHl), the physical pressure 0751), the density 
of the physical thermodynamic potential GB, the trace 
anomaly (1551) and the potential inhomogeneity (I5T1) as 
functions of volume V and temperature T for the free 
neutral scalar field on the finite lattice with the mass of 
7r° pion. For U = 3^ fm^ and fixed values of Np and 
Ncr, the function e(T)/T^ has two maxima, one maxi¬ 
mum at the lower temperatures and the other maximum 
at the higher temperatures. See the left panels of Fig. [5] 
This function vanishes as T —>■ oo. Thus, the lattice 
physical energy density (1751) does not recover the Stefan- 
Boltzmann energy density ess at high temperatures. 
The second maximum of the function e(T)/T'^ moves 
toward to higher temperatures with and the values 
of the function e{T)/T^ at the intermediate points be- 
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T (MeV) V (W) 


FIG. 2. (Color online) The exact lattice results. The physical energy density e (I78II . the physical pressure p (|79|l . the density 
of the physical thermodynamic potential ujph (|77l). the trace anomaly Aa dMI and the potential inhomogeneity Ai (HU as 
functions of volume V and temperature T for the free neutral scalar field with the mass of 7r° pion on the lattice for Njs = 20 
and No- = 20. Left panels: All curves were calculated at the volume V = 3^ fm®. Right panels: All lines correspond to the 
temperature T = 100 MeV. 


tween these two maxima tend to their Stefan-Boltzmann 
limit as N„- —>■ oo. The same behavior is clearly seen for 
the lattice physical thermodynamic potential ujph{T)/T'^. 
However, the lattice physical pressure p{T)/T‘^ has only 
one maximum at high temperatures and does not recover 
its Stefan-Boltzmann limit. The lattice trace anomaly 
A2(T) has a maximum at low temperature and vanishes 
with increasing T. However, the lattice potential inho¬ 
mogeneity —Ai(T) is an oscillating function of T and 
it is not equal to zero for those values of T for which 
the energy density and pressure are nonvanishing func¬ 
tions. Thus, as the lattice potential inhomogeneity Ai 
is not equal to zero, the statistical formalism for the free 
massive neutral scalar field on the lattice is thermody¬ 
namically inconsistent for all values of the temperature 
T at fixed values of A(g, and volume V. 

The volume dependence of the physical energy density 
d78l) , the physical pressure (17^ , the density of the phys¬ 
ical thermodynamic potential dzzl), the trace anomaly 
(l80l) and the potential inhomogeneity (IMI) is presented 
in the right panels of Fig. [71 The physical energy den¬ 
sity e{V) on the finite lattice increases infinitely with 
the volume V decrease. Such a behaviour is also seen 
in the density of the lattice physical thermodynamic po¬ 
tential ujphiy). However, the lattice physical pressure 


p{V) tends to zero as F —>■ 0. With growing volume V 
the function e{V) decreases, then increases and reaches 
a maximum, and drops to zero in the limit F —^ oo. 
The maximum of the function e(F) moves toward larger 
volumes F with increasing the lattice size and the 
values of the function e(F) at the intermediate points 
between the maxima of this function tend to a constant 
as iVo- —>■ oo. The same behavior is seen also for the 
lattice physical thermodynamic potential ujph{V). The 
lattice physical pressure p{V) increases with F, reaches 
a maximum and then drops to zero in the limit F —> oo. 
The lattice trace anomaly A2(F) and the lattice poten¬ 
tial inhomogeneity — Ai(F) increase infinitely as F —>■ 0. 
The trace anomaly A2 on the lattice is not equal to zero 
for all values of F for which the lattice physical energy 
density and the lattice physical pressure are nonvanishing 
functions. For large values of volume F > 100 fm^ the 
potential inhomogeneity Ai on the lattice is not equal to 
zero also for all values of F for which the lattice physical 
thermodynamic potential and the lattice physical pres¬ 
sure are nonvanishing functions. Thus, as the lattice po¬ 
tential inhomogeneity Ai is not equal to zero, the sta¬ 
tistical formalism for the free massive neutral scalar field 
on the finite lattice is thermodynamically inconsistent for 
all values of the volume F including the thermodynamic 
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FIG. 3. (Color online) Continuum limit. The physical energy density (1641) . the physical pressure (I68|l . the density of the 
physical thermodynamic potential (1711) . the trace anomaly A 2 (1721) . the potential inhomogeneity Aj (1731) and the mean 
number of particles (N)'^ as functions of volume V and temperature T for the free neutral scalar field with the mass of pion 
in the continuum limit. Left panels: the curves were calculated at the volume V = 3®, 6® fm® and in the thermodynamic limit 
(V — >■ 00 ). Right panels: the curves correspond to the temperature T = 100 and 200 MeV. 


limit. In the thermodynamic limit the physical pressure 
dZlD and the density of the physical thermodynamic po¬ 
tential dZZl) on the finite lattice vanish. 

Figure |3] represents the behavior of the physical en¬ 
ergy density (l64l) . the physical pressure (IM)l . the density 
of the physical thermodynamic potential (ED, the trace 
anomaly ( 1721 ) . the potential inhomogeneity ( 1731 ) and the 
mean number of particles (iV)'^ = ^ -1 / — 1) as func¬ 

tions of volume V and temperature T for the free neu¬ 
tral scalar field in the continuum limit with the mass of 
7r° pion. For 14 = 3^ fm^, the function e‘^{T)/T^ has 
a maximum at T ^ 30 — 40 MeV. With growing the 
temperature T this function decreases and attains a con¬ 
stant value as T —>■ 00 . For V = 6^ fm^, the function 


e‘^{T)/T'^ is a monotonically increasing function with T 
and tends to the same constant value as T —>■ c». The 
same behavior is seen for the physical thermodynamic 
potential —ujpf^{T)/T*; however, the maximum of this 
function is located at T ^ 50 MeV. For both V = 3^ fm^ 
and V = 6^ fm^, the physical pressure p^(T)/T^ increases 
monotonically with T and tends to a constant value as 
T —> c». The physical energy density e°, the physical 
thermodynamic potential and the physical pressure 
in the continuum limit tend to the Stefan-Boltzmann 
energy density esb, the Stefan-Boltzmann density of the 
thermodynamic potential ujsb and the Stefan-Boltzmann 
pressure psB, respectively, as T —>■ 00 . For both values 
of V, the trace anomaly A^iT) in the continuum limit 
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FIG. 4. (Color online) Continuum limit for the massless scalar field. The physical energy density ()64|l . the physical pressure 
p‘^ (1681) . the density of the physical thermodynamic potential (17111 . the trace anomaly A 2 (1721) . the potential inhomogeneity 
dial and the mean number of particles {N)‘^ as functions of volume V and temperature T for the free neutral scalar field 
with the mass m = 0 in the continuum limit without zero-mode term. Left panels: the curves were calculated at the volume 
1/ = 3® and 6® fm®. Right panels: the curves correspond to the temperature T = 100 and 200 MeV. 


has a maximum at T ^ 30 — 40 MeV. With growing the 
temperature T this function decreases and tends to zero 
as T -^ 00 . The appearance of the maximum in the trace 
anomaly in the continuum limit at low temperatures T is 
independent of the presence of the maximum in the en¬ 
ergy density. For both V = 3^ fm^ and V = 6^ fm^, the 
potential inhomogeneity — A)^(T) in the continuum limit 
has a maximum. The potential inhomogeneity — A);(T) 
for the free massive neutral scalar field in the continuum 
limit decreases with T and tends to zero as T —>■ 00 . The 
nonzero values of the potential inhomogeneity A); at fi¬ 
nite values of T indicate that for fixed values of V the 
statistical formalism in the continuum limit is thermody¬ 
namically self-consistent only in the limit T —>■ 00 . 

The volume dependence of the physical energy density 


(|M1). the physical pressure (1551) , the density of the phys¬ 
ical thermodynamic potential (El), the trace anomaly 
(El, the potential inhomogeneity (El and the mean 
number of particles {N^ in the continuum limit is pre¬ 
sented in the right panels of Fig. |3l The physical energy 
density in the continuum limit increases infinitely with 
the volume V decrease. Such a behaviour is also seen in 
the density of the physical thermodynamic potential in 
the continuum limit 0 ;°^. However, the physical pressure 
tends to zero as V —^ 0. With growing volume towards 
V > 10^ fm^ the physical energy density the physical 
thermodynamic potential and the physical pressure 
in the continuum limit reach their values in the ther¬ 
modynamic limit, which are different from the values of 
the Stefan-Boltzmann limit. The trace anomaly A 2 and 
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the potential inhomogeneity — AJ increase infinitely as 
y —>■ 0. The trace anomaly in the continuum limit 
is not equal to zero for all values of V and fixed temper¬ 
ature T. The potential inhomogeneity —for the free 
massive neutral scalar field in the continuum limit de¬ 
creases with V and vanishes in the thermodynamic limit 
as y —>■ oo. Thus, in the thermodynamic limit for large 
volumes and fixed values of T the statistical formalism 
for the continuum variables is thermodynamically self- 
consistent. 

Figure S] represents the behavior of the physical en¬ 
ergy density (l64l) . the physical pressure (l68l) . the density 
of the physical thermodynamic potential dnD, the trace 
anomaly ([72l) . the potential inhomogeneity (|73l) and the 
mean number of particles (A)'^ as functions of volume 
y and temperature T for the free massless neutral scalar 
field in the continuum limit without zero-mode term. For 
y = 3^ and y = 6^ fm^, the physical energy density 
e^{T)/T'^ coincides with the physical pressure p°(r)/T^ 
and they both are monotonically increasing functions 
with T. The same behavior is seen for the physical ther¬ 
modynamic potential however, the physical 

thermodynamic potential —u}pi^{T)/T‘^ does not coincide 
with the physical pressure p‘^(T)/T‘^. The physical en¬ 
ergy density £“, the physical thermodynamic potential 
ujpf^ and the physical pressure in the continuum limit 
tend to the Stefan-Boltzmann energy density the 

Stefan-Boltzmann density of the thermodynamic poten¬ 
tial UJSB and the Stefan-Boltzmann pressure pss, respec¬ 
tively, as T —^ oo. For all values of V, the trace anomaly 
A^iT) for the massless field in the continuum limit is 
equal to zero. For both y = 3^ fm^ and y = 6^ fm^, 
the potential inhomogeneity —A\{T) for the free mass¬ 
less neutral scalar field in the continuum limit has a max¬ 
imum. It decreases with T and tends to zero as T —> oo. 
The nonzero values of the potential inhomogeneity A^ at 
finite values of T indicate that for fixed values of V the 
statistical formalism in the continuum limit is thermody¬ 
namically self-consistent only in the limit T —>■ oo. 

The volume dependence of the physical energy density 
(|M1). the physical pressure (1551) , the density of the phys¬ 
ical thermodynamic potential (El), the trace anomaly 
dzl, the potential inhomogeneity dzi and the mean 
number of particles (A)'^ for the free massless neutral 
scalar field in the continuum limit is presented in the right 
panels of Fig.Hl The physical energy density coincides 
with the physical pressure and they both tend to zero 
as y —> 0. The physical thermodynamic potential in 
the continuum limit also tends to zero as y —>■ 0; how¬ 
ever, it does not coincide with the physical pressure 
With growing volume towards V > 10^ fm^ the physical 
energy density e°, the physical thermodynamic potential 
and the physical pressure in the continuum limit 
reach their values in the thermodynamic limit, which cor¬ 
respond to the values of the Stefan-Boltzmann limit. The 
trace anomaly A^ in the continuum limit is equal to zero 
for all values of V and fixed temperature T. The po¬ 
tential inhomogeneity — A^; for the free massive neutral 



FIG. 5. (Color online) The volume dependence of the mo¬ 
mentum gap. 

scalar field in the continuum limit decreases with V and 
vanishes in the thermodynamic limit as y —> oo. Thus, 
the nonzero values of the potential inhomogeneity A^ at 
finite values of V indicate that for fixed values of T the 
statistical formalism in the continuum limit is thermo¬ 
dynamically self-consistent only in the thermodynamic 
limit y —>■ oo. 

In the continuum limit the thermodynamic inconsis¬ 
tency (A^ ^ 0) of the system at low temperatures T and 
small values of volume V may be explained by the fact 
that in the volume V the system contains less than one 
particle on average, {NY < 1- See the lower panels of 
Figs. [3] and 0] Thus, the thermodynamic inconsistency 
of the system at low temperatures T and small values of 
volume y is an artifact of the grand canonical ensemble 
because in the grand canonical ensemble the term with 
the zero number of particles also contributes to the par¬ 
tition function. Note that for the free real scalar field 
there is no conserved current. However, the number of 
neutral particles can be accounted by the eigenvalues of 
the number operator 

N = J (Ypa'^{p)a{p). (93) 

The definition of this operator for the free real scalar 
field can be found, for example, in (^ . The mean num¬ 
ber of particles in the continuum limit {NY corresponds, 
namely, to the statistical average of this number operator 
A and has the meaning of the average number of neutral 
bosons in the volume V at temperature T [2^. With 
growing volume V and temperature T the mean number 
of particles in the continuum limit (A)'^ increases and the 
system described in the continuum limit becomes ther¬ 
modynamically self-consistent, i.e., A° —>■ 0. 

In a finite volume the momentum operator of parti¬ 
cles has the discrete eigenvalues. See Eqs. dH and (EH]). 
Not all values of momentum of particles are allowed. In 
the volume V the gap between the allowed values of the 
momentum of particles is Ap = Figure [S] repre¬ 

sents the volume dependence of the gap of momentum Ap 
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in the range of V characteristic for the ultrarelativistic 
heavy ion and hadron-hadron collisions. From Figs. [31 
m and [5] it is clearly seen that the free real scalar field 
is thermodynamically inconsistent (A^ ^ 0) for those 
values of V at which the momentum gap Ap is approxi¬ 
mately the same or larger than the temperature T of the 
system. We may conclude that the quantized momenta 
of particles in a small volume V and at finite tempera¬ 
ture T lead to the thermodynamic inconsistency of the 
free real scalar field in the grand canonical ensemble. 


VII. DISCUSSION AND CONCLUSIONS 

In the present paper, the grand canonical partition 
function (the thermodynamic potential) for the free neu¬ 
tral scalar field in one spatial dimension given in the 
framework of the method of path integral quantization 
was exactly derived in a general form by the new alge¬ 
braic method. The partition function in the momentum 
space was obtained from the partition function of the 
conhguration space by using the Fourier transform on a 
finite lattice. The matrix elements of the matrix A in 
one spatial dimension were found explicitly. The par¬ 
tition function in the momentum space was generalized 
to the three-dimensional spatial momentu m sp ace. It has 
recovered the partition function defined in (36l | . The main 
thermodynamic quantities of the grand canonical ensem¬ 
ble in three dimensions were obtained. In the continuum 
limit, the thermodynamic quantities for the free neutral 


scalar held given in the method of the path integral quan¬ 
tization resemble the second quantized quantities, which 
are obtained by the method of canonical quantization. 
We found that in the continuum limit the density of the 
thermodynamic potential of the free neutral scalar held is 
the sum of the vacuum and physical terms and does not 
contain any inhnite constant. The vacuum and physical 
terms of the thermodynamic quantities on a hnite lattice 
were dehned as in Ref. [3^ . 

The general condition for the thermodynamic self- con¬ 
sistency of any statistical system in the grand canonical 
ensemble was formulated. The thermodynamic proper¬ 
ties of the free neutral scalar held both on the hnite lat¬ 
tice and in the continuum limit were studied in the range 
of temperature and volume typical for the ultrarelativis¬ 
tic proton-proton and heavy-ion collisions. We found 
that the free massive neutral scalar held on the hnite 
lattice is thermodynamically inconsistent in the whole 
range of the variables of state T and V including the 
thermodynamic limit because its thermodynamic poten¬ 
tial is not a homogeneous function of the hrst order with 
respect to the extensive variable of state V. However, in 
the continuum limit the free neutral scalar held is ther¬ 
modynamically self-consistent only at high temperatures 
or in the thermodynamic limit. At hnite values of V it 
is thermodynamically inconsistent. 
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